Magnetic fields generated by r-modes in accreting millisecond pulsars 
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In rotating neutron stars the existence of the Coriolis force allows the presence of the so-called 
Rossby oscillations (r-modes) which are known to be unstable to emission of gravitational waves. 
Here, for the first time, we introduce the magnetic damping rate in the evolution equations of 
r-modes. We show that r-modes can generate very strong toroidal fields in the core of accreting 
millisecond pulsars by inducing differential rotation. We shortly discuss the instabilities of the 
generated magnetic field and its long time-scale evolution in order to clarify how the generated 
magnetic field can stabilize the star. 



PACS numbers: 04.40.Dg, 04.30.TV, 04.30.Db, 97.60. Jd 



INTRODUCTION 



R-mode oscillations are present in all rotating stars and 
they are unstable to emission of gravitational waves 
These modes play therefore a very important role in the 
astrophysics of compact stars and in the search for grav- 
itational waves. R-mode instabilities are associated with 
kinematical secular effects which generate differential ro- 
tation in the star and large scale mass drifts, particularly 
in the azimuthal direction. Differential rotation in turn 
can produce very strong toroidal magnetic fields inside 
the star and these fields damp instabilities converting the 
energy of the mode into magnetic energy. This mecha- 
nism has been proposed in the case of rapidly rotating, 
isolated and newly born neutron stars in Refs. [3,0]. In 
our work we consider the back-reaction of the magnetic 
fields on r-mode instabilities by inserting for the first time 
the magnetic damping rate into the evolution equations 
of r-modes. In this way we can follow the temporal evolu- 
tion of both magnetic fields and r-modes even on a long 
time scale. In particular we show that taking into ac- 
count the r-mode instabilities it is possible to generate 
very strong magnetic fields in accreting millisecond pul- 
sars. 

It is important to recall that r-mode instabilities are 
damped also by the viscosity of the system, either shear 
or bulk At low temperature a very important role 
could also be played by the so-called Ekman layer, lo- 
cated at the interface between the solid crust and the 
fluid inner core. The friction in the Ekman layer can be 
significantly enhanced with respect to friction in a purely 
fluid component. Anyway it is very difficult to give a pre- 
cise estimate of this effect [H [f| , and in our analysis we 
have not included the Ekman layer in order to study in a 
more clear way the effects of magnetic fields on r-modes. 
We will anyway shortly discuss how the presence of the 
Ekman layer can modify our results. 
The paper is organized as follows: in Sec. II we write the 
r-mode equations taking into account also the magnetic 



damping rate; in Sec. Ill we discuss the equations regu- 
lating the time evolution of the magnetic field and of the 
damping rate; in Sec. IV we discuss numerical solutions 
of the equations; in Sec. V we discuss open problems 
concerning the possible presence of the Ekman layer, the 
effect of a superconducting shell and instabilities of the 
generated internal magnetic field. Finally in Sec. VI we 
drive our conclusions. 



II. R-MODE EQUATIONS 

We want to derive the equations regulating the evolu- 
tion of r-modes in the presence of a pre-existent poloidal 
magnetic field B p and of the generated internal field. We 
start by considering the conservation of angular momen- 
tum, following Ref. [6]. The total angular momentum J 
of a star can be decomposed into an equilibrium angu- 
lar momentum J* and a canonical angular momentum J c 
proportional to the r-mode perturbation: 



j = j*(m, n) + (i-Kj)j c 



J c = -K c o?J, (1) 



where Kij tC \ are dimcnsionless constants and J* = 7*f2. 
Following Ref. [7[ the canonical angular momentum 
obeys the following equation: 



dJ c /dt 



2J c {F g (M,n) 
-[F v (M,il,T v ) 



■F mi (M,Sl,B)]} (2) 



where F g is the gravitational radiation growth rate, F v — 
F s + Fb is the sum of the shear and bulk viscous damping 
rate, T v (t) is a spatially averaged temperature and now 
we also have introduced the damping rate F mi , associated 
with the generated internal magnetic field. The explicit 
expression of F mi will be discussed in the following. 
The total angular momentum satisfies the equation: 



dJ/dt = 2J c F g + j a (t) - hQF n 



(3) 



'Electronic address: cuofano@fc.infn.it 



where J a is the rate of accretion of angular momen- 
tum (we have assumed it to be J a — M(GMR) 1 ^ 2 , see 
Ref. (fl) and F ms is the magnetic breaking rate associated 
to the external poloidal magnetic field. In the present pa- 
per we are not considering further breaking mechanisms 
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III. MAGNETIC DAMPING 
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FIG. 1: Top panel: instability region in the Temperature 
vs Frequency plane obtained using either the bulk viscos- 
ity damping rate Ft given in Jl| (region above the dot- 
dashed line) or the Fi, given in [S| (region above the dashed 
line). Also shown are the paths followed by accreting neu- 
tron stars, without toroidal magnetic fields (solid lines). The 
three curves correspond to different values of accretion rate 
M = (10 -8 , 10" 9 , 10" 10 ) Mb yr" 1 . We plot also the tem- 
perature equilibrium curve [3] (dotted line), obtained taking 
into account the neutrino cooling and the reheating due to 
viscosity. Bottom panels: new paths obtained taking into ac- 
count the new generated toroidal fields. Here Bd = 10 8 G. 
The moments at which toroidal magnetic fields damp r-mode 
instabilities are indicated by asterisks. 



as e.g. the interaction between the magnetic field and 
the accretion disk. Combining Eqs. ([2]) and ([3]) we ob- 
tain the evolution equations of the r-mode amplitude a 
and of the angular velocity of the star Q: 



da 
~dt 



a{F g ~F v -F mi ) 

+a[K 3 F g + (1 - K 3 )(F V + F mi )]K c a 2 



aM ( G 



1/2 



aF„ 



dt 



2in \MR 3 J 2 
-2K c tla 2 [K 3 F g + (1 - Kj){F v + F m J] 
MQ M ( G , ' " 

1 lF 7r 



M I \MR 3 



(4) 



(5) 



where = IMR 2 with / = 0.261 for an n=l polytrope 
and K c — 9.4 x 10 -2 , see Ref. @. Our results turn out to 
be rather insensitive to the value of Kj ~ 1 (see Ref. 0). 



The crucial ingredient introduced before is the mag- 
netic damping rate F mi , which we have inserted in the 
evolution of r-modes. The expression of the magnetic 
damping rate has been derived in 0, H|, where it has 
been shown that while the star remains in the instabil- 
ity region, the r-modes generate a differential rotation 
which can greatly amplify a pre-existent magnetic field. 
More specifically, if a poloidal magnetic field was origi- 
nally present, a strong toroidal field is generated inside 
the star. The energy of the modes is therefore transferred 
to the magnetic field and the instability is damped. 
We assume that the stellar magnetic field B is initially 
dipolar and aligned with the star's spin axis 



B = B*(i - 0) 



R 3 

Bd—^(2 cos 9 e r 



sinOee) (6) 



where Bd is the strength of the equatorial magnetic field 
at the stellar surface. The initial magnetic field is as- 
sumed to be described by Eq. ^ both inside and out- 
side. In order to avoid divergences we assume the previ- 
ous radial dependence of the dipolar field to hold only for 
r > B/2. If for instance the dipolar field remains roughly 
constant in the inner region its contribution to the total 
energy of the final toroidal field is negligible. Therefore 
in the following and in agreement with Ref. [3|], radial 
integrations extend from R/2 to R. 

To estimate the magnetic field produced by r-modes we 
start by writing the I = m = 2 contribution to the per- 
turbation velocity: 



6v(r,6,<t>,t)=anR{-) Yge M . 



(7) 



Following Ref. we get the total azimuthal displace- 
ment from the onset of the oscillation at to up to time t, 
which reads: 



AS^i) ee / Sv^it^dt' 

= l(^)k 2 (9)j\ 2 (t')n(t / )dt / + 0(a 3 ) 



(8) 



where k 2 {9) = (l/2) 7 (5!/7r)(sin 2 6» - 2cos 2 6>). The re- 
lation between the new and the original magnetic field 
inside the star in the Lagrangian approach reads Qj: 



B j B k dx-(t) 

— (x,t) = — (x,t ) fc 3 ;/ 

p p dx k (t Q ) 



(9) 



This equation implies that the radial dependence of the 
initial and final magnetic field is the same. Integrating 
on time the induction equation in the Eulerian approach 
one gets 0]: 



8B° 



6B r ~Q 

B e [ 4>{t')dt' ~ B e [ S -^ldt' (10) 
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where is the toroidal component. 

The expression of the magnetic damping rate reads: 



(dE M /dt) 
E 

4(1 -p) 



BpRA' J Q * a 2 (t')Q,(t')dt' 



9irp- (8.2 x to- 3 ) 



Mfl 



(11) 



where E is the energy of the mode, Em is the magnetic 
energy, A' » Oil) is a dimensionless parameter [3[ and 
p = 0.5 indicates that we are integrating from R/2 to 
R. The time integral over the r-mode amplitude a takes 
contribution from the period during which the star is 
inside the instability region. 



IV. NUMERICAL SOLUTIONS 

In our numerical analysis we use the estimate given in 
[l[ for the gravitational radiation reaction rate due to the 
I = m = 2 current multipole 



F g = ^M lA R\ Q Pli s" 1 



(12) 



as well as for the bulk and shear viscosity damping rates 
1 



F b = 
F s = 



^-^MdRnPlil* s- 
1 



2.2 x 10 7 



Mi.4-R 10 5 T 9 2 s 1 



(13) 



where we have used the notation Mi. 4 = M/1.4 Mq, 
R 10 = R/10 Km, P_ 3 = P/l ms and T 9 = T/10 9 K. 
It is important to remark that the damping rate Fb de- 
pends on the density profile of the star. Using a different 
profile in Q they obtained a value for FJ, three orders of 
magnitude larger than in [l[ . Fortunately our results are 
insensitive to this quantity, only the instability window 
is modified as shown in Fig. [TJ 

Viscosity depends critically on temperature. We in- 
clude three factors in modeling the temperature evolu- 
tion: modified Urea cooling, shear viscosity reheating, 
and accretion heating. The cooling rate due to the mod- 
ified Urea reactions, e u , reads [l(| 



e u = 7.5 x 10 39 Ml'^Tl 



erg s 



(14) 



The neutron star will be heated by the action of shear 
viscosity on the r-mode oscillations. The heating rate 
due to shear viscosity, e s , reads [l[ 



= 2o?VL 2 MR 2 JF t 
= 8.3 x 10 3 Vft 2 JM*£R W 



j a. f 9/4 p- 15/4^-2 



erg s 



(15) 



where J = 1.635 x 1CT 2 . 

Accretion heating has two components. We use the es- 
timates given in [ll| . The first contribution arises when 
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FIG. 2: Top panel: temporal evolution of r-modes amplitude 
without toroidal fields. Middle panel: temporal evolution of 
r-modes amplitude with toroidal fields for Bd — 10 8 G (solid 
line) and Bd = 10 G (dotted line). Bottom panel: temporal 
evolution of volume averaged toroidal magnetic fields. We 
consider M = (10~ 8 , 10~ 9 , 10 -10 ) M yr _1 . 



accreted matter undergoes nuclear burning at the surface 
of the star 

e„ = X 1.5 MeV = 4 x 10 51 Mi 4 erg s _1 (16) 

m B 

where ms is the mass of a baryon. 

The second contribution arises because the flow is as- 
sumed to be advection dominated. The heating rate is 

^-T^ = 8xl ° 5lM i 3 4 /6 ^i.4 ergs" 1 . (17) 

Finally we use the estimate of the heat capacity C v given 
in [U 

C v = 1.6 x 10 39 M^T g erg K _1 . (18) 
The equation of thermal balance of the star is 
1 



d 



:C V T 



(19) 



We consider a scenario in which the mass accretion spins 
an initially slowly rotating neutron star up to millisec- 
ond period and we investigate the evolution of internal 
toroidal magnetic fields when the star enters the r-modes 
instability window. 

We start by discussing the evolution of temperature and 
spin frequency obtained without magnetic fields. In 
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Fig. Q] we show that the star crosses the r-mode insta- 
bility region in a regime of thermo-gravitational runaway 
Q. R- modes grow exponentially due to the decrease of 
the shear viscosity with increasing temperature. As a 
consequence r-mode amplitudes rapidly reach the satu- 
ration value (we chose a sa t = 1 although, as we shall see 
later, magnetic fields limit a <C a sat ) and the viscosity 
heats significantly the star. At this stage the star loses 
angular momentum by emission of gravitational waves 
and goes out of the instability region in a time of hun- 
dred of years. 

Taking into account magnetic fields, the evolutionary sce- 
nario for the star is quite different. We discuss results ob- 
tained solving in a self-consistent way Eqs. (|4I5|11|19|) . 
From a practical viewpoint we proceed as follows: we 
first solve Eqs. (|4I5I19|) and we get an estimate of F mi (t) 
inserting the results in Eq. (| 1 1 p : we use this estimate to 
solve again Eqs. (|4l5ll9p . This procedure is iterated till 
numerical convergence is reached. 

The exponential growth of r-modes induces extremely 
large secular effects and the toroidal magnetic field is 
either produced or amplified by the wrapping of the 
poloidal field produced by the secular velocity field which 
is mostly toroidal. In the bottom panel of Fig. [T] we 
show the new trajectory of the star in the Temperature- 
Frequency plane obtained taking into account the gen- 
erated toroidal fields, and we indicate with asterisks the 
moments at which magnetic fields damp r-modes. It is 
important to remark that this happens when the star is 
still in the region which was unstable taking into account 
only the viscous damping. 

In Fig. [2] we show the evolution of the r-mode ampli- 
tude a, without magnetic field (top panel) and with mag- 
netic field (middle panel). Three different values of ac- 
cretion rate M and two values of the initial poloidal mag- 
netic field Bd are considered. In the scenario with mag- 
netic fields, the maximum values of a are in the range 
Umax ~ [10" 6 — IO" 1 ] and the generated toroidal fields 
are in the range B tor ~ [10 11 - 10 15 ] G (Fig. H bottom 
panel). In this analysis we have considered values of ac- 
cretion rate M and magnetic field Bd typical of accreting 
Low Mass X-Ray binaries (LMXBs). 
In Fig. [3] we show the generated toroidal magnetic field 
as a function of the temperature at which the star enters 
the instability window. We display the volume averaged 
field obtained at the moment in which the field itself has 
completely damped r-modes. The value of the magnetic 
field is the asymptotic one displayed in the bottom panel 
ofFig.H Results ofFig.[3]are obtained using Eqs. I|16ll7p 
to describe the connection between accretion rate M and 
temperature T, but they are rather insensitive to the ex- 
act relation between these two quantities. 



V. OPEN PROBLEMS 

In our work we have not discussed the possible exis- 
tence of Ekman layer. If this layer is present it stabilizes 
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FIG. 3: Volume averaged toroidal magnetic field generated 
by r-modes as a function of the temperature at which the 
star enters the instability window. For reference we show also 
the value of the temperatures which correspond to a mass ac- 
cretion rate M = (IO -8 , IO" 9 , IO" 10 ) M yr" 1 . The various 
curves correspond to different values of the initial poloidal 
magnetic field, in the range [10 s — 10 9 ] G. 



the star up to frequencies of a few hundred Hertz. The 
scenario discussed above does not change qualitatively, 
but since the star enters the instability region at higher 
frequencies, the growth rate of magnetic fields is larger. 
Another important open point concerns the possible for- 
mation of superconductivity in the core of neutron stars. 
A fraction of the core, whose temperature is below the 
critical value T c ~ 10 9 K, is expected to be either a 
Type I QJ, [l| or a Type II superconductor [3, Q]|. The 
exact nature of the possible superconducting layer is still 
uncertain as well as the precise value of the supercon- 
ducting gap. The fundamental quantity in our analysis 
is the size of the superconducting region which we show 
in Fig. fusing the results obtained in [l4|. The thickness 
of the superconducting shell is about 1 km. In the same 
figure we also display the energy density of the r-modes. 
It is clear that most of the volume where r-modes can 
develop remains not superconducting and in that region 
our analysis can still be applied: r-modes are damped by 
the magnetic field and therefore the perturbation velocity 
drops to zero. R-modes can develop in the thin supercon- 
ducting layer, but their amplitude is strongly suppressed 
because the perturbation velocity has to vanish at the in- 
terface between the superconducting and the not super- 
conducting region. Moreover, if the superconductivity 
is of Type II, the magnetic flux of the original poloidal 
field organizes into quantized flux tubes. R-modes will 
stretch these flux tubes, increasing their length and thus 
changing their magnetic energy. As in the case of the nor- 
mal core, this process can generate very intense magnetic 
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FIG. 4: Energy density of r-modes as a function of the star 
radius. The estimates about the size of the superconducting 
region are quite uncertain. Recent results obtained in (141 ] 
suggest that the superconducting shell has a thickness < 1 
km and it is located in the grey area. 



fields which lead to the damping of the r-mode oscilla- 
tions even in the superconducting layer Q ■ 

Several issues remain open concerning how the new 
generated magnetic fields are affected by possible insta- 
bilities. In the stably stratified environment of a stellar 
interior there are two types of instabilities: the Parker 
(or magnetic buoyancy) and the Tayler instabilities (or 
pinch- type), both driven by the magnetic field energy in 
the toroidal field. The buoyancy instability is negligible 
for B tor J$ 10 15 G, so we focus on Tayler instabilities 
because they set in at a lower field strength [l6|. It is 
important to remark that in a stable and stratified neu- 
tron star the condition for the Tayler instability reads 
0: 



UJA_ 

n 



> 



9. 



1/2 



V 
r 2 il 



1/4 



(20) 



where loa = B /(47rp) 1 / 2 r is the Alfven frequency, ~ 
5 x 10 4 s -1 is the compositional contribution to the buoy- 
ancy frequency and r\ ~ 1CP 9 cm 2 s _1 is the magnetic 
diffusivity |18|. From Eq. (|2T)|) we can conclude that in 
the stably stratified core of a neutron star, the Tayler 
instability sets in for BZ r > 10 12 G. 
After the development of the Tayler instability, the 
toroidal component of the field produces, as a result of 
its decay, a new poloidal component which can then be 
wound up itself, closing the dynamo loop. Both compo- 
nents then grow, more slowly, until the saturation level is 
reached, when the field is being destroyed by the instabil- 
ity at the same rate at which it is being amplified by the 
differential rotation fl9| . When the differential rotation 
stops the field can evolve into a stable configuration of 
a mixed poloidal-toroidal twisted-torus shape embedded 
inside the star with an approximately dipolar field con- 
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FIG. 5: New r-modes instability regions for neutron stars with 
stable configurations of mixed poloidal-toroidal fields in the 
inner core (no magnetic field - solid line; Bt or , P oi = 10 12 G - 
dashed line; Btor, P oi = 10 13 G - dotted-dashed line; Btor, P oi = 
10 14 G - dotted line). We assume that poloidal and toroidal 
components have similar strengths. 



nected to it outside the star [20|, lU [22|, [23| • Once the 
field is stabilized it should evolve as a result of diffusive 
processes as Ohmic dissipation, ambipolar diffusion, and 
Hall drift [24[, whose typical time-scales are: 



tohmic 
tambip 
tflall 



3 x 10 M 



5 x 10 s 



vLl ( 
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n\ 


Pnuc 




y r 






p 


B 12 \ 


Pnuc 



yr 



yr 



(21) 
(22) 
(23) 



where L§ = L/10 cm is the size of the region embedding 
the magnetic field. 



VI. CONCLUSIONS 

We have shown how r-modes can generate strong 
toroidal fields in the core of accreting millisecond neutron 
stars, and how these fields influence the growth rate of r- 
mode instabilities. Tayler instability sets in for strengths 
of the generated fields of the order of 10 12 G and stabi- 
lizes the toroidal component by producing a new poloidal 
field of similar strength. This stable configuration evolves 
on a time-scale, regulated by diffusive processes. Our re- 
sults imply that in the core of accreting neutron stars in 
LMXBs, rotating at frequencies v > 200 Hz, there are 
strong magnetic fields with strengths B > 10 12 G. 
Finally, it is tempting to try to investigate how the new 
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stable configuration of magnetic fields modifies the insta- 
bility window of r-modes. More explicitly the scenario we 
have in mind is the following: 

• the compact star inside the LMXB initially enters 
the instability window and it follows the trajecto- 
ries that are displayed in the bottom panel of Fig. [I] 
and in Fig. O When the toroidal field reaches the 
critical value dictated by the Tayler instability the 
poloidal field suddenly increases and we assume it 
reaches a value of the same order of B^ r . The new 
magnetic configuration, in which the poloidal and 
the toroidal field are of the same order of magni- 
tude, is stable. In the new configuration the in- 
ternal poloidal field is a few orders of magnitude 
larger than the initial poloidal field. R-modes try 
again to deform the poloidal field generating a new 
toroidal component but the magnetic damping rate 
given by Eq. (fTTj) is now much larger and the star 
is stable with respect to r-modes up to frequencies 
of the order of 200 Hz. 

• the star continues accreting and increasing its fre- 
quency and it enters again the instability region 
generating again a new toroidal field. When the 
new toroidal component becomes much larger than 
the poloidal component Tayler instability sets in 



again increasing the value of the poloidal field. In 
Fig. [5]we show how the instability window changes 
in dependence of the internal magnetic configura- 
tion, assuming that the toroidal and poloidal com- 
ponents are equal. 

• to determine the actual maximum value of the lim- 
iting frequency we take into account also the time- 
scale associated with the diffusive processes de- 
scribed by Eqs. ([2"l"l[2"3]) . If the magnetic fields ex- 
ceed a few 10 14 G the diffusion time-scale becomes 
shorter than 10 6 yr and the star cannot accelerate 
further even for very large values of M. 

It is clear that even in the absence of the Ekman layer 
the new internal magnetic fields can stabilize stars with 
frequencies up to several hundreds Hertz. The possible 
presence of a superconducting shell would screen these 
internal fields so that they would not affect the dynamics 
of the external region. The phenomenological implica- 
tions of this model will be discussed in a forthcoming 
paper. 
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Pagliara and Andreas Reisenegger for many useful and 
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